Successive opening of the Fermi surface in doped A^-leg Hubbard ladders 
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We study the effect of doping away from half-filling in weakly (but finitely) interacting A'^-leg 
Hubbard ladders using renormalization group and bosonization techniques. For a small on-site 
repulsion U , the A^-leg Hubbard ladders are equivalent to a N-haxiA model, where at half-filling the 
Fermi velocities are vi = vn < V2 = vn-i < . . .. We then obtain a hierarchy of energy-scales, 
where the band pairs (j, N + 1 ~ j) are successively frozen out. The low-energy Hamiltonian is then 
the sum of N/2 (or (A'^ — l)/2 for A'^ odd) two-leg ladder Hamiltonians without gapless excitations 
(plus a single chain for A'^ odd with one gapless spin mode) — similar to the A'^-leg Heisenberg 
spin-ladders. The energy-scales lead to a hierarchy of gaps. Upon doping away from half-filling, the 
holes enter first the band(s) with the smallest gap: For odd A^, the holes enter first the nonbonding 
band {N + l)/2 and the phase is a Luttinger liquid, while for even A'', the holes enter first the band 
pair {N/2, N/2 + 1) and the phase is a Luther-Emery liquid, similar to numerical treatments of the 
t-J model, i.e., at and close to half-filling, the phases of the Hubbard ladders for small and large U 
are the same. For increasing doping, hole-pairs subsequently enter at critical dopings the other band 
pairs (j, N + 1 — j) (accompanied by a diverging compressibility): The Fermi surface is successively 
opened by doping, starting near the wave vector (7r/2,7r/2). Explicit calculations are given for the 
cases A" = 3, 4. 

PACS number(s): 71.10.Pm 
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I. INTRODUCTION 

Interest in the field of ladder systems grew rapidly, 
when the possibility of superconductivity in doped two- 
leg ladder materials was proposed [0,^. Experimentally, 
such a behavior has been observed in a two-leg ladder 
material under high pressure |^] . The A'^-leg ladders have 
then been studied both as a first step towards the two- 
dimensional (2D) case Q and because of their unusual 
physical properties. For example, spin-ladders exhibit 
an odd-even effect [|] ; while even-leg ladders have a spin- 
gap, odd-leg ladders have one gapless spin mode. Sim- 
ilarly, the behavior upon doping is very different. The 
two-leg ladder is a Luther-Emery liquid, but the lightly 
doped three-leg ladder consists of a Luttinger liquid (LL) 
plus an insulating spin liquid (ISL) |^,^. The latter is 
due to a truncation of the Fermi surface (FS) in two chan- 
nels. Such a truncation is argued to take place also in 
the 2D case, where upon doping some parts of the FS 
become conducting and other parts remain an ISL |^. 

Analytic works for (doped) ladders are usually based 
on the Hubbard model and restricted to a weak on-site 
repulsion U (for large U, see Ref. |9|). Insulating behav- 
ior at half-filling is then due to umklapp processes. The 
two-leg ladder has been investigated by many authors 
over the past few years [p^-p^. We note that in the 
lightly doped two-leg ladder paired holes form a dilute 
gas of hard-core bosons, which are equivalent to spinless 
fermions ||ll|,|l^ . Away from half-filling, the phases of the 
weak interaction limit [/ ^ of the three- and four-leg 
ladder have been derived in Refs. (for possible su- 

perconductivity, see Ref. ||l6|). For the three- leg ladder, 



a C2S1 phase has been obtained upon doping (a phase 
with n gapless charge- and m gapless spin modes is de- 
noted by CnSm). This contrasts numerical findings for 
the (strongly interacting) t-J model, where at low dop- 
ings the phase is ClSl [^j^. An idea to explain this 
discrepancy has been given in a related context, where a 
one-dimensional electron gas (IDEG) was coupled to an 
insulating "environment" (another IDEG) [|l7|l . 

In this paper, we study the doping S away from half- 
filling of A''-leg Hubbard ladders in the case of weak but 
finite repulsive interactions, < U <^ tj_,t {t and tj_ 
denote the hopping matrix elements along- and between 
the chains). Previous works have taken the limit U ^ 
first, keeping the doping S finite Here, we treat 

the opposite limit, i.e., we keep U finite and then inves- 
tigate the lightly doped case S ^ 0, where the umklapp 
interactions present at half-filling cannot be neglected. It 
is then advantageous to take the half-filled A^-leg ladders 
as a starting point. 

For a small U, the Hubbard-model reduces to a weakly 
interacting A^-band model. After linearization around 
the Fermi points, each band is characterized by a Fermi 
velocity vj; at half-filling vi = vn < V2 = wjv-i < . . .. 
The low-energy physics is then derived by a combination 
of the renormalization group (RG) and bosoniza- 
tion methods [|l9|,^. Integrating the RG equations 
(RGEs) at half-filling, we find a decoupling into band 
pairs [j, N + 1 — j) and a hierarchy of energy-scales 
Tj ^ te~°''"^/^ {a is a constant of the order of 1), where 
the band pairs scale towards a two-leg ladder fixed point 
and become frozen out (this behavior was not noted in 
previous works in particular, Ref. concen- 
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trated on doping levels (well) away from half- filling). For 
N even, all excitations are then gapped. For odd N, the 
remaining band behaves like a single chain such that the 
charge degrees of freedom acquire a gap while the spin- 
excitations stay gapless. We therefore recover the odd- 
even effect that is present in Heisenberg spin-ladders. 

The gaps for the band pairs are given by Aj ^ Tj re- 
sulting in a hierarchy of gaps. Upon doping, the holes 
enter first the band(s) with the smallest gap. For odd 
N, the gap of the single chain is the smallest one and at 
lowest doping, the phase is a LL (ClSl). For even N, 
the spin degrees of freedom remain gapped upon doping 
and the phase is a Luther-Emery hquid (CISO). These 
phases are the same as obtained in (numerical) treat- 
ments of the strongly interacting t-J model [||J|,0,^,^, 
i.e., at and close to half-filling, the phases for small and 
large U are the same. For increasing doping, paired 
holes enter at critical dopings, where each time the com- 
pressibility diverges, subsequently the other band pairs 
{j,N -I- 1 — j). In other words, the FS is successively 
opened, first near the wave vector (7r/2,7r/2), spreading 
towards the wave vectors (tt, 0) and (0, tt) for increasing 
doping, see Figs. |l| and ^ Explicit calculations will be 
given for the cases = 3, 4. 

In Sec. II, the A^-leg Hubbard model is introduced 
and the phases at half-filling are derived. The main re- 
sult of this section is a low-energy Hamiltonian allowing 
us to study the effect of doping away from half-filling: 
In Sec. Ill, including the doping as a perturbation of the 
half-filled Hamiltonian, we treat first the case when only 
one band (pair) is doped and then the case when many 
band pairs become doped. 



II. HALF FILLED TV-LEG HUBBARD LADDERS 

A. The Hamiltonian 

The A^-leg Hubbard model is given by H = Ha + Hi^t , 
where the kinetic energy is 

HQ = ~tJ2 4s + ^)djsix) + H.c. 

-ix^4+i.(a;)rf,s(a;)+H.c. (1) 

and the interaction term is 

Hint = uY,d]^{x)d,^{x)d]^{x)d,y{x). (2) 

The hopping matrix elements along- and perpendicular 
to the legs are denoted by t and t± respectively and U > Q 
is the on-site repulsion. For weak interactions, U <ti t±,t, 
it is advantageous first to diagonalize Hq, 



Ho=J2J dke,{k)^l{k)'^js{k), (3) 

where ^t^j^, and ^'^^ are the creation- and annihilation 
operators for the band j. 
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FIG. 1. The bands of the three-leg ladder correspond to 
three different points on a 2D FS (denoted by 1, 2, and 
3). The square is the umklapp surface, which is the FS at 
half-filling (for t = t±). The half-circles indicate the opening 
of the FS, when the band two is doped. At (close to) the wave 
vector (■7r/2, 7r/2), the phase is then a LL, while at (37r/4, 7r/4) 
and (7r/4, 37r/4), the phase is an ISL. 

For open boundary conditions perpendicular to the legs, 
the dispersion relations ej are given by (j = 1, . . . , N) 

ej{k) = -2i cos(fc) - 2t±_ cos [nj/iN + 1)] . (4) 

The Fermi momenta in each band kpj are determined 
by the chemical potential fi — £j{kpj) and the filling 
n = 2(7riV)^^ '^kpj- Since we are only interested in the 
low-energy physics, we linearize the dispersion relation 
around the Fermi momenta, resulting in Fermi veloci- 
ties Vj ~ 2tsm{kFj). We introduce operators '^ji/Ljs for 
right- and left movers at the Fermi level in each band. 
At half-filling n = 1, /i = and the velocities are given 

by 

Vj ^vj^2 \lt^-{ti_ cos[7rj7(iV + l)]}^ (5) 
where j — N + \ — j . Note that 

Vi=VN <V2^ VN-l < ■■■■ (6) 

These particular values of the Fermi velocities will lead 
to a hierarchy of energy scales (see below). 

For generic Fermi momenta, the interactions are 
forward- (/) and Cooper- (c) scattering within a band 
or between two different bands. Following Ref. the 
interacting part of the Hamiltonian can then be written 
as 



2 




FIG. 2. The bands of the four-leg ladder correspond to 
four different points on a 2D FS (denoted by 1, 2, 3, and 
4). The square is the umklapp surface, which is the FS at 
half- filling (for t = t±). The half-circles indicate the open- 
ing of the FS, when the bands two and three are doped, 
forming a Luther-Emery (LE) liquid, where the pairing takes 
place between quasi-particles at (close to) the wave vec- 
tors (27r/5,37r/5) and (Stt/S, 27r/5). At {An/5,n/5) and 
(tt/S, 47r/5), the phase is still an ISL. Similarly as for the 
two-leg ladder, the order-parameter has a different sign in the 
Fermi points 2 and 3. 



-ffint — ^ J dx {f[jJRiiJLjj — fij^Rii ■ J 



■ ''Ljj 



+ ^ j dx {c'^jjJRjjJLjj - CJjJrjj ■ 3ljj) , (7) 

where the U(l) and SU(2) currents are defined as 

Jh^, ^Y.'^ls'^hjs (8) 



and 



k 

his^ss' ^hjs' 1 



(9) 



where again h = R/L and the r*^ are the PauU matri- 
ces (fc = x,y,z). Due to symmetry, f^-" = j''f and 

In the half-filled case we must generalize the interact- 
ing Hamiltonian (^. The Fermi momenta of the bands 
J and J add up to tt, kpj + kp^ — tt, allowing umklapp 
interactions to take place between these two bands (for 
"jJJj differ in notation from Ref. by a factor 2), 

(4j^IJ + H.C.) - (4^- . + H.C.) 

+"^^;I;(4.^^^ + 4l^..+H-c.). (10) 
We have used the definitions 



and 



(11) 



(12) 



where e = — ir^. For odd A^, there exists a single-band 
umklapp term for the band r = (A^ + l)/2 



[iRrr^Lrr + H.C. 



(13) 



In addition to these processes, there exist various inter- 
actions involving four different bands or three different 
bands (only for N odd). For the half-filled three-leg lad- 
der, three different three-band umklapp processes 

"1223 (42^23 + 432^21 + H.C. 
+"2213 (422^13 + ^13^22 + H.C. 
-"?223 (lm2 • Ii23 + 1^2 ' 1^21 + H.C.) (14) 

and two different three-band non-umklapp processes 

^1223 {Jr12Jl23 + Jr32Jl21 + H.C.) 
-Ci223 (Ji?12 • Jl23 + ^ R32 ' Jl21 + H.C.) (15) 

take place and have to be added to Eq. (^. 

In the four-leg case, the four-band umklapp interac- 
tions read 



•'1234 



4l2-^i34 + Irs^Ilu + H.C. + ^^324 (2 ^ 3) 



1423 



(^m4-^i23 + lR23^Lli + H. 



(I!, 



-"1234 1 ■l_R12 ■ 1^34 + ljj34 ' Il12 + H.C, 



X324 (2 --3) (16) 



and the four-band non-umklapp interactions take the 
form 

Cl234 {JRl2JL3i + Jra?,Jl21 + H.C.) + C1324 (2 ^ 3) 
"^1234 (Jfll2 • Jl34 + Jfl:43 ' Jl21 + H.C.) 

-C?324 (2 3) . (17) 



In our case vi = V4 and V2 = W3 implying C4234 = c5'324, 
"1234 = "1324 ■ The generalization to TV > 4 is straight- 
forward such that we do not reproduce it here. 

The (initial) values of the couplings are of the order of 
the on-site repulsion U (for N = 3, see Appendix A). 
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B. Low-energy phases 

Using the RG method (and bosonization), we derive 
the low-energy phases at half-filhng and a Hamiltonian 
allowing us to study the effect of doping, i.e., doping is 
then included as a (small) perturbation. We have gener- 
alized the results given in Refs. ||ll|,|l^ in order to treat 
ladders with N > 2 legs at half- filling For weak 

interactions, ?7 <C it is sufficient to take the RGBs in 
one-loop order, i.e., quadratic in the coupling constants. 

Integrating the RGEs, we find, that the Fermi ve- 
locities Vj [see Eqs. (||) and (0)] lead to a hierarchy of 
energy-scales by Tj ~ te~°'^^'^, where the band pairs 
(jiJ) become successively frozen out. This exponential 
dependence on vj/U is a consequence of the logarithmic 
derivatives of the couplings with respect to the energy in 
the one-loop RGEs. As we show below, the low-energy 
Hamiltonian is then the sum of N/2 [{N - l)/2 for N 
odd] two-leg ladder Hamiltonians corresponding to the 
band pairs (j, j) (plus the Hamiltonian of a single chain 
for TV odd). 



(21) 



where, using Eq. (|18|), we have expressed the interactions 
of band 1 and 3 in terms of (713. The eigenvalues of M 
are (7, —1, —1, —1) and 



f gisindl' ^ ^In 
vi + V3 Ji^ 12 



513(0 



513(^0) 



(22) 



where Iq < I denotes the scale at which the interactions 
of the bands 1 and 3 are sufficiently close to the two-leg 
ladder ratios. The three-band interactions therefore stay 
in fixed ratios and 
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(23) 



such that the three-band interactions are arbitrary small 
at the scale where 313 ~ t. 



1. Three-leg ladder 

It is instructive to start with the RG flow of the three- 
leg ladder in the limit vi — v:} V2 (i.e., t±/t ^ V2), 
where we can neglect in the RGEs all contributions, 
where contractions over the band 2 take place (see Ap- 
pendix A). We then obtain the following results. First, 
the renormalization of the two- and single-band interac- 
tions of the bands 1 and 3 is exactly the same as for a 
two-leg ladder, where the couplings flow towards univer- 
sal ratios (for a comparison, see Ref. ITT 
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0, 



(18) 



(19) 



and become of the order of the bandwidth t at the tem- 
perature scale T\ ~ te""^"^/'^ (ai is a constant of the 
order of 1). Second, the interactions between the bands 1 
and 2 (and 2 and 3) are either not renormalized or flow to 
0. Finally, the scaling of the remaining three-band inter- 
actions can be calculated as follows (note that U2213 stays 



always small), 
and 



M = 



Defining v = (q223'C?223' 



^1223' "1223 
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(20) 



we rewrite the RGEs as (the scaling variable I is related 
to the energy by T ~ te~^^ , see Appendix A) 




logU/t 



FIG. 3. The figure shows (for t = t\_) the strong cou- 
pling value as a function of the initial value U /t, where we 
fix the strong coupling value of C13 at the bandwidth t. The 
three-band coupling C1223 then decreases to (the flow of the 
other three-band couplings is similar), as the initial value is re- 
duced, and the coupling approaches its asymptotic value, 
cjg = C13/4. We note that also for the half-filled two- leg lad- 
der [11], the flow to the asymptotic ratios is very slow, i.e., 
one has to take very small initial values U /t. 

For other ratios of the velocities vi/v2 < 1, we have 
performed a numerical integration of the RGEs: Plotting 
the strong coupling value as a function of the initial value 
U/t (we fix the strong coupling value of cf^ at the band- 
width t), we find that the couplings given in Eq. ( |T^ ) 
always grow up to the bandwidth t, approaching their 
universal ratios, while all the other couplings stay an or- 
der of magnitude smaller, i.e., the asymptotic behavior 
of the couplings for [//t — > is the same as for vi ^ V2 
(see Fig. |^). We note that in Ref. |T^ a flow to strong 
coupling of the three-band interactions has been found 
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for t±/t < 0.86; we only obtain that for vi w V2, i.e., in 
the strongly anisotropic limit t±/t < 0.2 pi] . 

To conclude our RG analysis of the three-leg ladder: 
The couplings of the bands 1 and 3 scale towards the two- 
leg ladder ratios, become of the order of the bandwidth t 
at the temperature scale Ti ~ te~°'^'"^^^ , and at (and be- 
low) Ti, the bands 1 and 3 are decoupled from the band 2. 
The Hamiltonian becomes therefore H — H13 + Hg. The 
first term H13 is the Hamiltonian of a two-leg ladder (see 
Ref. and the second term is the Hamiltonian of a 
single "chain". At the scale Ti, the couphngs of H13 are 
of the order of t and all the charge and spin-modes of the 
bands 1 and 3 acquire a gap (this result is obtained by 
bosonization, see Ref. The couplings of Hg are still 

small at Ti, i.e., < 222,1*22 ^ t, such that the charge 
and spin- modes of band 2 remain gapless. 

For a (small but) finite U, we then have to investigate 
the system at energies below Ti . For the Hamiltonian Hg 
it is well-known that it has a charge gap below a temper- 
ature T2 ~ te-°'2^2/u (since a2 « ai we write in the fol- 
lowing a for both of them) and that the spin-excitations 
remain gapless [Q. Since the two- leg part has no gap- 
less mode, the final phase is COSl — the same as for the 
strongly interacting hmit, i.e., a three-leg spin-ladder. 

2. Four-leg ladder 

For the four-leg ladder, we find a similar behavior. The 
couplings of the bands 1 and 4 grow, become of the order 
of the bandwidth t at the energy scale Ti ~ te~°''"^^^ , 
and stay in the same ratio as for a two- leg ladder, 
while all the other couplings remain small. In the limit 
vi = V4 V2 — V3, the RGEs of the four-band couplings 
become the same as for the three-band couplings and con- 
sequently, similarly as above, the four-band interactions 
stay in fixed ratios and 

<23J<, « iU/4,f'' ■ (24) 

At the temperature scale Ti , the bands 1 and 4 are there- 
fore decoupled from the bands 2 and 3. Again, for a fi- 
nite U, we then have to study the flow of the couplings 
of band 2 and 3 at energies below Ti . We find that they 
flow at the scale T2 ~ te'"""^^^ also to a two-leg lad- 
der fixed point (the flow to the two-leg ladder fixed point 
takes place for a wide range of initial values, also for non- 
Hubbard-like). The (bosonized) Hamiltonian is then the 
sum of two two-leg ladder Hamiltonians, H = Hi4 + H23- 
All charge and spin excitations are therefore gapped (as 
it is the case in the Heisenberg four-leg spin ladder). 

3. The case N > 4 

Analyzing the RGEs for iV > 4, we conjecture that 
this successive freezing out of band pairs (plus a 



single band for N odd) at the characteristic energies 
- ie"""^/'^ takes place also for (smaU) iV > 4 with 
U ^ t/N. Again, the Tj are a result of the one- loop 
RGEs: the couplings of the band pair scale towards 
a two-leg ladder fixed point and become of the order of 
the bandwidth t at Tj . Note that we have a hierarchy of 
energy-scales 

Ti>T2> ...>%, (25) 

where for N even r = N/2 and for N odd r = {N + l)/2. 

To conclude this section: At low energies the Hubbard- 
ladder Hamiltonian becomes for N even the sum of N /2 
two-leg ladder Hamiltonians and for TV odd the sum of 
{N — l)/2 two-leg ladder Hamiltonians plus the Hamil- 
tonian of a single chain 

H = ^ Hjj + (5jv,oddffs, (26) 
j 

where for N odd, the single chain [the band (A^-fl)/2] has 
the lowest energy-scale and for even the two-leg ladder 
Hamiltonian Hjq/2^N/2+i [corresponding to the band pair 
{N/2, N/2 + 1)]. The (half-filled) two-leg ladder Hamil- 
tonians have no gapless excitations and the single-chain 
Hamiltonian present for odd N has only one gapless spin- 
mode. The phases at half-filling are thus the same as 
for the Heisenberg spin-ladders [|j. Since the Hubbard 
model converges for large U onto the t-J model (which 
is at half-filling equivalent to the Heisenberg model) [^ , 
the phases of the half-filled A^-leg Hubbard ladders are 
the same for small and large U. 

We would like to emphasize that we have obtained the 
final result, Eq. (p6|), describing decoupled band pairs, 
from an explicit analysis of the coupled RGEs of the in- 
teracting A^-band problem. 

III. HOLE DOPED LADDERS 

Since the low-energy cutoff is typically of the order of 
te"*/'^, the effect of the umklapp interactions has to be 
taken into account for low dopings, S < e"*/*^. It is then 
advantageous to take the half-filled ladder as a starting 
point to study the effect of doping. Similarly as done in 
Refs. [ p"2p^ ], we model the doping by adding a chemical 
potential term —fiQ, where Q is the total charge. 

In a single chain, a hole is decomposed in two quasi- 
particles, one with charge but no spin (holon, here the 
fermionic "kink" ) and the other with spin but no charge 
(spinon). The quasi-particles of a two-leg ladder are 
bound hole-pairs (singlet of charge 2). At low doping, 
these hole-pairs behave as hard-core bosons, which are 
equivalent to spinless fermions. In Sec. HI A, we show 
that the lightly doped A^-leg ladders can indeed be de- 
scribed by an effective A^/2-band [{N -h l)/2-band for A^ 
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odd] model of spinless fermions. We then study the case 
when only one of these bands is doped. In Sec. Ill B, we 
increase the hole doping such that many bands become 
subsequently doped. 



A. One band (pair) doped 

For the following discussion, we use the hosonized form 
of the Hamiltonian (|2^). Bosonization is the 

method of rewriting fermionic creation (annihilation) op- 
erators in terms of field operators and 11^ satisfying 
the commutation relation a{x) ,Ha{y)] = i^ix — y) (a 
labels the charge and spin modes for the different bands; 
for the application of bosonization to ladders in the limit 
U <C t_L,i, see, e.g., Refs. pl| , p^ , p6| ) . It is convenient to 
introduce the dual field of dxOa — Hq. For fermions 
on a chain or ladder, bosonization applies in the contin- 
uum limit. 

It is straightforward to generalize the effective Hamil- 
tonian for the half-filled two-leg ladder obtained in 
Ref. to N/2 [{N - l)/2 for N odd] two-leg ladder 
Hamiltonians. We introduce the fields 



$pj) and n+j = Bp 



(27) 



combining the charge fields of the band pair (j, j) . Drop- 
ping bare energy terms, the Hamiltonian then reads 



g+j cos(v87r$+j) - 2fj,dx^ 



(28) 



plus for odd N the Hamiltonian of a single chain (drop- 
ping the gapless spin-part) 



bands by Aj 
Eq. (H)) 



Tj, resulting in a hierarchy of gaps (see 



Ai > A2 > . . . > 



(31) 



The equations of motion following from the Hamilto- 
nians (^^ and (^) are a set of sine-Gordon (SG) equa- 
tions. Here, the physical solutions of the SG equations 
are the soliton-solutions. Solitons are fermionic particles 
(kinks). The kinks of the single-chain field $ps have a 
mass Ms A^ and represent the charge part of single 
holes, while the kinks of the two-leg ladder fields <f>+j 
have a mass Mj ~ Aj and represent paired holes (sin- 
glets) . The band structure of these solitons reads 



(32) 



At low doping 5 < e^*/^, the half-filled A^-leg ladders are 
therefore described by an effective iV/2-band [(-/V-l- l)/2- 
band] model of spinless fermions, where the bottom of 
the bands are at the energy-scales Aj (for iV = 3, see 
Fig. I). 
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— (a,$p,)^ + i^p,H^, 



-gps cos(\/87r$ps) - fid^'^ps^- 



(29) 



The coupling g^^-j is of the order of gjj, the Ka are the 
Luttinger liquid parameters (at half-filling, K^j 1 and 
Kps 1/2 1 13 2^), and the Ua are the velocities of the 



charge modes. The charge (doping) is given by 



(30) 



We note that (^) and ( |2^ ) are the usual Hamiltoni- 
ans describing a commensurate-incommensurate transi- 
tion The fields of the band pair (j,j), (and of 
the single chain $ps), become pinned at the energy-scales 
Tj ~ te~"''"'/^ . We therefore define effective gaps for the 



FIG. 4. The lightly doped three-leg ladder can be described 
by a two-band model of spinless fermions, where the lower 
band II contains the charge part of holes and the upper band I 
paired holes (singlets). At half-filling the phase is COSl; upon 
doping, the charge part of the holes enters the band II and the 
phase is ClSl. For dopings 5 > 5c (here, nSc = \kc\), paired 
holes enter also the band I and the phase becomes C2S1. 

At half-filling, /i = and minimizing the energy gives 
$+j K, K, 0. When doping the first hole (pair), the 
chemical potential /i = /i((5) jumps from to A^ and a 
soliton-solution with charge Q > minimizes the energy, 
i.e., the lowest- lying effective band of spinless fermions 
becomes doped. For odd N , the lowest-lying band is the 
single-chain band {N -I- l)/2 and the phase upon dop- 
ing is a dilute gas of spinous and kinks (ClSl), while 
for even iV, the effective band A^/2, corresponding to the 
two-leg-ladder band-pair {N/2, N/2 + 1), is doped first 
and the phase consists of singlets of paired holes (CISO). 
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Similarly as for the two-leg ladder, the (singlet) pair-field 
operator 



(X exp (^—iV^TrOpj^ cos(\/27r$crj) 
has a different sign in band j and j pif, 



(33) 



(34) 



such that the CISO phase has a cZ-wave like symmetry. 

It is instructive to rewrite the band operators 'i^hjs re- 
spectively the pair-field operators Pj in terms of the chain 
operators dhjs- For = 3, 



^h2s = "^(^''1^ ~ dh3s) 



and we find that the holes are situated on the outer 
legs. Both the phase and the location of the holes is 
in agreement with previous numerical treatments [^,0. 
For = 4, using P2p| ~ —1, we obtain 

P2 cx 0.22 (rffli|(iL2i + dR2]dLii + R ^ L) 
+0.22 {dR3idL4i + dBA'\dL3i + R<-^ L) 
-0.36 idRi^dL4i + diii^dLii + R ^ L) 
-0.14 {dR2^dL3l + dR3^dL2i +R^ L), (36) 

such that the singlets are on the top two legs, the bottom 
two legs, on the legs 1 and 4, and with lowest probability 
on the legs 2 and 3, similarly as found in Ref. Q for the 
t-J model. Finally, also the phases for N = 5,6 are in 
agreement with numerical works jj]. Since the t-J model 



is the large U limit of the Hubbard model 1 25 1 , we con- 
clude that at and close to half-filling the phases of the 
iV-leg Hubbard ladders are the same for small and large 
U. For the two-leg ladder, such a "universal" behavior 
has already been noted in Ref. |13|. 



B. Many bands doped 

When increasing doping, the chemical potential /i — 
li{S) increases too and the hole-pairs enter subse- 
quently also the higher-lying (effective) bands of spinless 
fermions. We thus obtain a series of critical dopings 5cj , 
where the first hole-pair enters the band j. For = 3, 4 
Sci is estimated as 



Since for S = 6cj, dfi/dS = 0, the compressibility 

(38) 



d6_ 
9/1 



diverges. The phase transitions at 6cj belong to the same 
universality class as the commensurate-incommensurate 



transition [g^. We note that our qualitative estimate 
for Scj does not allow for a comparison with numerical 
works. 

The effective band j corresponds to the band pair 
(j,j), such that for a given doping, the band pairs 
(A^/2, N/2 -f 1), . . . , are conducting, while the band 
pairs (j — 1, J-t- 1), . . . , (1, N) form still an ISL. Next, we 
interpret this result by mapping each band (pair) on a 
2D FS, see Figs. || and ||. Using the dispersion relation 
(^), the longitudinal Fermi momentum of the band j is 
given by 



kpj = TT — arccos 



■ cos 



N + 1 



(39) 



(35) and the corresponding transverse Fermi momentum reads 



kpyj — 



N+1 



(40) 



Upon doping, the holes enter (for t± ^ t) then first near 
the wave vector (7r/2, 7r/2) and finally for increasing dop- 
ing also near (tt, 0) and (0, tt). The FS thus becomes 
truncated, similarly as in the 2D case studied in Ref. [|. 

For increasing doping, a crossover takes place to the 
situation away from half-filling, where the umklapp pro- 
cesses can be neglected. It is instructive to study the 
phases for N = 3, 4, when both effective bands of spin- 
less fermions are doped and to compare with the phase 
away from half-filling. In the case iV = 3 (see Fig. |[), the 
different charges in the effective bands I and II forbid 
scattering processes of the form 



*flj*i?//*i/*L// + 1 ^ II, 



(41) 



where ^ n/j^i, creates a spinless fermion in band b = I /II, 
since they break U(l) invariance and only the interac- 
tions 



I ^ II 



(42) 



are allowed. For weak interactions, we can bosonize Eq. 
(^ ) resulting in a density term of the form dx^idx^ii 
(plus the same for the dual fields 9ijf). This term does 
not reduce the number of gapless modes such that the 
phase for (5 > (5c is C2S1. The same phase has been 
found by the RG and bosonization methods away from 
half- filling ||l^,|l^. This contrasts numerical calculations 
for the strongly interacting t-J model 0, where a spin- 
gapped phase has been found away from half-filling. Such 
a phase maybe appears also for weak interactions, if the 
condensation to the low-energy phase takes place for all 
bands at the same energy-scale (for another idea, see 
Ref. 0). 

For = 4, we have two types of paired holes with 
the same charge in the ladder, localized along the rungs. 
In contrast to the three-leg case, scattering processes of 
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the form (|l]) are allowed. We have shown in Ref. [E8|, 
that in such a two-band model of spinless fermions with 
repulsive interactions, slightly away from the band edge, 
binding takes place and the C2S0 phase close to the band 
edge becomes a 4-hole CISO phase. Away from half- 
fiUing, by KG and bosonization, the phase C3S2 has been 
found . Similarly as already noted in Ref. ||T^ , we ar- 
gue that the C3S2 phase is a result of the U ^ limit 
considered in Ref. ; and indeed, our phases agree with 
a numerical treatment of the t-J four- leg ladder At 
low doping, the authors obtained a phase of paired holes 
and at higher doping levels, four-hole clusters (or two- 
pair clusters). We therefore argue that for increasing 
doping the phase remains [/-independent in even-leg lad- 
ders, but not in odd-leg ladders. 

We finally note that including next nearest neighbor 
hopping terms etc., the Fermi momenta of band pairs add 
up no more exactly to tt, i.e., Vj ~ Vj. However, since the 
phases are controlled by the relative size of gaps, only 
a sufficiently large perturbation can qualitatively change 
the corresponding physics. 



(effective) band j corresponds to the band pair (j, j) 
and the bottoms of these bands are at the energy-scales 
Aj . Upon doping, the hole-pairs then subsequently enter 
these bands and the FS is successively opened, first near 
the wave vector (7r/2, 7r/2) and finally also near (7r,0) 
and (0,7r), similarly as in the 2D case Q. In addition, 
we obtain peaks in the compressibility at critical dopings 
dcj, where the first hole-pair enters the band pair (j,j). 

For increasing doping, the [/-independence of the 
phase seems to break down, at least for odd-leg ladders. 
However, the [/-independence and the similarity to the 
2D case at and close to half-filling support the hope that 
also in 2D, some of the weak-coupling RG results can be 
extended to large U S. 
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IV. CONCLUSIONS 

We have studied the doping away from half-filling in 
TV-leg Hubbard ladders in the case of a small on-site re- 
pulsion U , where the ladders are equivalent to a iV-band 
model. At half-filling, the N Fermi velocities are such 
that vi = vn < V2 ^ vn-1 < ■ ■ ■■ Using RG and 
bosonization techniques, we have obtained a hierarchy 
of energy-scales Aj ~ te~°''"'/^ , where the couplings of 
the band pair (j,j) (and of the single band {N + l)/2 
for N odd) are of the order of the bandwidth t, decou- 
pled from the other bands, and at a two-leg ladder fixed 
point. The low-energy Hamiltonian is then the sum of 
N/2 [{N -l)/2 for N odd] two-leg ladder Hamiltonians, 
corresponding to the band pairs (j,j), without gapless 
excitations (plus for odd N the Hamiltonian of a single 
chain with one gapless spin- mode). We thus obtain the 
same phases as for the Heisenberg spin-ladders. 

Upon doping, the holes enter first the band (pair) 
with the smallest gap Aj. In odd- leg ladders, this 
is the nonbonding band {N -\- l)/2 and the phase is 
a dilute gas of holons (here, the fermionic kinks) and 
spinous (ClSl), while in even-leg ladders, the band pair 
{N/2, N/2 + 1) is doped first and the phase consists of 
singlets of paired holes (CISO), reflecting the odd-even ef- 
fect at half-filling (i.e., of Heisenberg spin-ladders). The 
same phases have been found in numerical works for the 
t-J model [g,||j^j2l[l, i.e., at and close to half-filhng, the 
phases of the iV-leg Hubbard ladders are the same for 
small and large U. 

For a finite U and at low doping, the ladders can be 
described by an effective Af/2-band model [(TV -I- l)/2- 
band model for TV odd] of spinless fermions, where the 



APPENDIX A: RENORMALIZATION GROUP 

The renormalization group (RG) method is a con- 
trolled way of subsequently eliminating (integrating out) 
high-energy modes in a given Hamiltonian. While the 
noninteracting part Hq is at a i?G fixed point, couplings 
of the interacting part Hint may grow or decrease un- 
der a RG transformation (i.e., when lowering the en- 
ergy) . A subsequent (perturbative) elimination of high- 
energy modes in Hint then results in RGEs, which give 
the change of the coupling when lowering the energy, see 
Ref. [0. 

Below, we give the RGEs for the couplings of the three- 
leg ladder. Fully taking into account the symmetries [i.e., 
SU(2), U(l), and vi = v^], there are 21 independent 
RGEs. Using operator product expansion for the prod- 
ucts of the various currents, the derivation is straight- 
forward l|2^. Similar RGEs are obtained in the case 
TV > 3 

For the half- filled three- leg ladder vi — v^, such that 



P:0- _ rp,(J p,cr _ 
J23 ' ^ 



C23 : and c1\_ 



^33 



The initial 



values of the two and single band couplings are given by 



3U 



U 



U 



16 



'12 ~ ^12 



' ^22 



J13 ^ '-II ~ ''IS ~ 1 c ' •' : 



^1331 



,P _ 
"! '^1133 



32 



; ""1331 — ""22 ^ 



(Al) 



and fZ = 4 ff^ , cf, = 4cf,- . The initial values of the three- 

" ''J " ''J ''J ''J 

band couplings read 



rP - 
-1223 — 



u 



u 



u 



u 



'■1223 



"j "2213 ~ Q I "1223 



(A2) 
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The RGEs take the following form, where the energy- 
scale is related to ^ by T ~ te~'^K Starting with the 
above initial values the RGEs are then numerically inte- 
grated up to the scale where the first coupling is of the 
order of the bandwidth t. For one-loop RGEs, Ic = at/U, 
where a ^ 1 (this can easily be checked by multiplying 
the RGEs with t/U'^). We use the abbreviations 



A+ = ± 



('-1223 ) + 



^ (C?223)^' 



+ (""1223) + (^1223) 

= i4 (c5'223Cl223 ^ '^1223^1223) 
~ (^1223) ~ (""1223) 



-13 



p,(T 



2ui 



fP,0 n P 

J 13 _|_ ^^12 



1^12 



+ 



2V2 



D = 



^1223^1223 ~ fi- 



le 



1223 "1223 > 



(A3) 



where V12 = Vi + V2- 

The two-band forward-scattering interactions are 
renormalized according 



df?3 


1 


dl 


~ 2^ 
1 


<2 


1 


dl 


V\2 




1 


dl 


~ 2t^ 
1 


4ff2 


1 


dl 


fl2 



K3r + ^K3)' + K33l)' 
16 «133)^ + 4 K33l)^ 



(c?2)' + ^(cr2)'+4(«f2i3r 



16 

2 



+ 



2V2 



+ 



V\2 



2C13C13 — - (C13) — (/Ji 



2^1331^1331 



1 / cr x2 

2 1^1331^ 



2^12^12 2 ^^12) (/12) 



the single-band interactions according 



A-V2 

2V12' 



(A4) 



dl 



+ 



k^l 
1 

2v^ 



2vk 



3 

16 



dl 



^ 2^;fe 



(^1331) + ("1331) 



fc5^2 



+ 1 (.^2)^ + ^ 

■^2 t": 



= -E- 



fe#l 
1 

'2u^ 



2vk 



(cifc) + 2c^j.Cij. 



2t;i 



(c?i)' 



2'Wl33 1^1331 + 2 



cr -,2 
1331 J 



22 



dl 



E- 



(C2fe) + 2c^fcC2i; 



the two-band back-scattering interactions according 

+ ~ ^^13/13 + YQ'^lsfia + 4Uii33Mi33i~j + 



dc{2 
dl 



dl 



dc^2 
dl 



2 

V12 



+ 7 ( Ci2/f2 + ^2*^12/12 ) + 7r~'"2213^1223 



4 

1'12 



+ ~ ^^13/13 + fiis/fs 2^^13/13 + 4u^j33Uj33i^ 



1 



— E 9,„ ( '^lfe'^fc2 + ClftCfe2 + o'^lfe'^fc2 



( C12/12 + fil2/f2 ~ o ''12/1^2 
''12 V ^ 

4 „ , 



V12 



^2213 "1223' 



(A6) 



and for the three-band interactions, we find 

dCi223 _ '"l33l'"l223 + 4tt]'i33M]'223 ^^^33]^ 71^223 



dl 



2vi 



32wi 



97;'' 7;'' '^/-'^ 

^"22 "1223 , /^P 1 '-"-1223 /-fcr 
T Ci223<-^- "I ~ — 



dc' 



1223 



t^2 



^1223 p 



16 



H r <-i223 

'1331 + 4u5'l33 ~ o'"l331 



2vi 



+■ 



*1331"-1223 



2vi 



2^22^1223 
V2 



'-1223 

+Cl223^- + ^1223^- 2~^+" ^''^'^^ 

The two- and single-band umklapp interactions are renor- 
malized according 

^^1331 _ "1331 (fP ,„P\, ^"'1331 / ra _ „o- \ 
^Z~~^^^13+CllJ + ^^(/l3 

_^f£l3<133 ^ D 
Vl V2 

du-iiss _ c^aUi^^i + 4/f3tt ^^33 ^ ScJgU^^gg]^ ^ D 



2vi 



32vi 2v2 



1331 _ "1331 /fa 



dl 



(/f3-cf,) + ^(/f3 + cf,) 

■"1331 / f<7 I „CT \ I 4Ci3uJ^33 
Ul3 + CllJ 



2wi 



Vl 



1 



+ 2^2 ^''1223'"l223 + 201223^^1223 2cf223'*^1223) 



dU22 2022^22 _|_ ^ 

dl V2 Vl 



(A8) 
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and the three-band umklapp interactions according 



""1223 
dl 



n.P cP + 4?;'' c'' 
"l33l"^1223 ^ ^"ll33''1223 



3u 



2vi 



1331'-^1223 
32^1 



4u: 



2213'"12 
W12 



2u?, 



22'"1223 
V2 



'^"2213 
dl 

'^^1223 

dl 



I „p 'J"l223 ^cr 

+"1223"-"+ 

^''^1223'^12 I 3^^223*^12 



1^12 



'"1223 

2wi 



8W12 



^"2213-'12 
W12 



'1331 



4u?i33 



'■1331 



U 



-U 



1331'"1223 
1223^- ^ 



*"2213'^12 



2w; 



22'-1223 



'1223 



■1'12 



^^2 
^1223 r^G 



(A9) 
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